Effects of Random Link Removal on the Photonic Band Gaps of Honeycomb Networks 
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We explore the effects of random link removal on the photonic band gaps of honeycomb networks. 
Missing or incomplete links are expected to be common in practical realizations of this class of 
connected network structures due to unavoidable flaws in the fabrication process. We focus on the 
collapse of the photonic band gap due to the defects induced by the link removal. We show that the 
photonic band gap is quite robust against this type of random decimation and survives even when 
almost 58% of the network links are removed. 



Photonic crystals, and in particular photonic band gap 
(PBG) materials, dil constitute a class of dielectric ma- 
terials, in which the basic electromagnetic interaction is 
controllably altered over certain frequencies and length 
scales. [H, Hf The ability to tailor the photonic density of 
states in a prescribed manner and to engineer the symme- 
try properties of the electromagnetic field inside a PBG 
material enables the design of materials and devices that 
control the flow and the emission and absorption proper- 
ties of electromagnetic radiation. d-U Most of the pho- 
tonic crystals fabricated to date possess intrinsic disor- 
der, including size variation among the scattering centers, 
roughness of the dielectric walls, point defects, lattice 
dislocation and stacking defaults. foL fioj The presence of 
disorder results in defects modes that can fill the PBG 
region partially or completely and alter significantly the 
photonic properties of the dielectric crystals. 

In this paper we focus on the effect of random link re- 
moval on the properties of photonic gaps of honeycomb 
networks. It is well established that two-dimensional di- 
electric structures based on connected network architec- 
tures favor the formation of large PBGs for transverse 
electric (TE) polarized radiation. fi"ll - [il| The largest 
known TE PBG is achieved in a dielectric structure ob- 
tained by decorating the links of a honeycomb network 
lattice with dielectric walls of a certain thickness [l4| 
(also depicted in the inset of Fig. [TJ. 

We consider a special type of disorder produced by the 
random removal of links from a honeycomb connected 
network. This type of disorder is expected to be com- 
mon in practical realizations of this class of connected 
network structures due to unavoidable flaws in the fab- 
rication process. The major goals of this study are to 
explore the robustness of the PBG against random re- 
moval of links and to establish a relationship between 
the topology of defect complexes (resulting from the re- 
moval of multiple links) and the size of the band gap. 



The class of photonic structures considered in this 
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FIG. 1: (Color online) Band structure for the optimal unper- 
turbed TE PBG structure consisting of dielectric walls with 
dielectric constant e = 11.56 thickness w/a = 0.94 (shown as 
an inset), which presents a TE PBG of Aoj/ljc = 52.63%, the 
largest known PBG in two-dimensional photonic crystals for 
this dielectric constant contrast ratio. 

study are generated by decorating the sides of a two- 
dimensional honeycomb tessellation with dielectric walls 
of dielectric constant e = 11.56 (corresponding to Si in 
the infrared region) and thickness w/a = 0.94, where a 
is the honeycomb lattice constant, corresponding to the 
largest known TE PBG: 52.63% (see Fig.[TJ). 

To evaluate the influence of the link disorder on the 
band gap properties, we construct a finite domain of the 
honeycomb dielectric network and define a statistically 
independent probability, p, for the removal of the dielec- 
tric links (i.e., Bernoulli percolation [lj| [IH) that con- 
nect the vertices of the honeycomb network. For each p, 
a statistical ensemble consisting of 1500 configurations is 
generated, and the photonic band structure is calculated 
along a closed contour around the Brillouin zone (as in 
Figs. 1, 3a and 4a) for each member of the ensemble. 
In addition, we have checked that gap size by calculating 
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FIG. 2: (Color online) The small filled circles represent the 
PBG size as a function of the percentage of randomly removed 
links, as obtained by averaging over a statistical ensemble of 
1500 configurations. The upper and lower lines represent the 
maximum and minimum value among the ensemble. The inset 
shows the extrapolation of the average PBG to percentages 
large enough to close the photonic band gap. 



the band structure for 100 k-points randomly distributed 
in the first Brillouin zone. Due to computational limita- 
tions, the domain considered consists of 8 x 8 unit cells. In 
Fig. [2 we plot the maximum, minimum and average PBG 
in a given ensemble of randomly generated structures as 
a function of the number of links removed. We note that 
while full connectivity is associated with the largest PBG 
structures [l7j , the band gap is very resilient against line 
removal; even after the random removal of p = 16.7% of 
all the links, the structures present a sizable average band 
gap of Au/wc = 17% (here Aw is the PBG size and ojc 
the gap central frequency). Moreover, the minimal PBG 
configuration still retains a gap of about 8%, while the 
maximal PBG configurations display sizable band gaps 
of about 28%. 

We have analyzed different types of defects generated 
by the link removal and their influence on the band gap. 
In general, the link removal results in a lowering of the lo- 
cal dielectric constant, which, as expected, generates an 
increase in the frequency of the photonic modes around 
the lower edge of the PBG. These modes ascend inside 
the PBG and subsequently the PBG size is reduced. As 
the density of the removed links increases, the unoccu- 
pied links may cluster together to generate larger defect 
complexes. Isolated clusters containing a small number 
of missing links generate cavity-like defects. Cavity-like 
defects induce photonic modes with a high degree of lo- 
calization and low dispersion, and generate defect bands 
with small spectral bandwidth. 

In Fig. [3] a we plot the band structure for the configu- 
ration that displays the largest gap (Aw/wc = 27.78%) 
in a structure with p = 16.7%. We note that the modes 
shift to inside the PBG due to the presence of material 




FIG. 3: (Color online) The connected network photonic struc- 
ture with the largest TE PBG for p = 16.7% (selected from 
1500 randomly generated configurations), a) The photonic 
band structure for dielectric constant contrast e = 11.56 
and wall thickness w/a = 0.94, displaying a PBG with 
Au/wc = 27.78% (shaded region), b) The magnetic field 
distribution for the mode at the lower edge of the PBG (the 
black lines show the contour of dielectric material distribu- 
tion). 



defects present minimal dispersion. The corresponding 
dielectric structure (the black-line contour in Fig. [3] b) 
and the magnetic field distribution of the photonic mode 
defining the lower edge of the PBG confirm that most of 
the defect modes are confined around the random cavity- 
like defects in the structure. 

A second type of defect is a waveguide-like complex, 
which occurs when removed links concatenate to create a 
defect-cluster spanning the entire structure. In this case, 
the defect modes are associated with propagation along 
the waveguide channel, thus localized only transverse to 
the waveguide. As a result, they display stronger disper- 
sion and have a more deleterious effect on the width of 
the PBG. 

In Fig. 2] a, we plot the band structure for the configu- 
ration that achieves the minimal PBG (Aut/uic = 7.96%) 
among structures with p = 16.7%. Here we note that 
among the modes that partially fill the PBG region, there 
are modes showing considerable dispersion. As shown in 
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FIG. 4: (Color online) The connected network photonic struc- 
ture presenting the smallest TE PBG for p = 16.7% (se- 
lected from 1500 randomly generated configurations), a) 
The photonic band structure for dielectric constant contrast 
e = 11.56 and wall thickness w/a = 0.94, displaying a PBG 
Aw/wc = 7.95%. b) The magnetic field distribution for the 
mode corresponding to the lower edge of the PBG (the black 
lines show the contour of dielectric material distribution). 



Fig.[4]b, the unoccupied links form a waveguide channel 
(the dielectric profile is shown by the black- line contour), 
and the mode at the lower edge of the PBG has the mag- 
netic field rather tightly confined inside the waveguide. 

To summarize, we have analyzed the effects of link dis- 
order on the TE band gaps in honeycomb connected net- 
work photonic structures. We have demonstrated that 
the PBG is robust against this type of disorder and have 
identified the relevant defect topologies. Removing links 
of dielectric material introduces zero-dimensional (cav- 
ities) or one-dimensional (waveguide channels) defects, 
and, depending on their spectral location and dispersion, 
the PBG size is reduced as p increases. At some threshold 
value of p, the PBG vanishes. Thus, our study can be re- 
garded to be a photonic percolation problem. We define 
the photonic percolation threshold as the value p = p c 
at which the PBG first vanishes in the infinite-system 
limit. For p > p Cl the PBG is identically zero. Compu- 
tational requirements make difficult a precise determina- 
tion of the photonic percolation threshold, but we have 
performed an extrapolation of our results (shown in the 
inset of Fig. [2]), which indicates that the PBG vanishes if, 
on average, p c = 57.6% of the links are removed. Further 
band structure calculations confirm that, if 57.6% of the 
connected network links are randomly removed, the aver- 
age PBG collapses. This photonic-percolation-threshold 
value should be contrasted with that of Bernoulli nearest- 
neighbor bond percolation on a honeycomb network for 
which p c = 1 - 2sin7r/18 w 65.28%. The fact that the 
values of p c are not identical is not surprising because 
they depend in detail on different physical characteristics. 
Other percolation phenomena with yet different values p c 
include elasticity percolation [iH lla | and quantum per- 
colation. [l|| 
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